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Abstract. All the observations available in the Chandra and XMM-Newton archives have been used to investigate the X-ray 
spectral properties of 3C 33. In this paper is presented a complete X-ray analysis of the nuclear emission of this narrow line 
radio galaxy. The broad band spectrum of 3C 33 is complex. The hard part resembles that of Seyfert 2 galaxies, with a heavily 
obscured nuclear continuum (N H ~10 23 cirT 2 ) and a prominent Fe Ka line. This represents the nuclear radiation directly 
observed in transmission through a cold circumnuclear gas. On the other hand an unabsorbed continuum plus emission lines 
seem to fit well the soft part of the spectrum (0.5-2 keV) suggesting that the jet does not significantly contribute to the X-ray 
emission. We discuss the possible collisional or photoionized origin of the gas that emits the soft X-ray lines. Our results, 
strengthened by optical spectroscopy favor the photoionization scenario. 
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1. Introduction 

With the advent of high resolution grating spectrometers, 
X-ray spectroscopy has become the most powerful tool for 
investigating the ionized component of the gaseous AGN 
environment. The detection of emission lines in Seyfert 2 
galaxies has revealed gas photoionized by the nuclear engine 
(Sako et al. 2000; Kinkhabwala et al. 2002; Bianchi et al. 
2005; Pounds & Vaughan 2006; Guainazzi & Bianchi 2007). 
Spatially resolved studies realized with the Chandra telescope 
on at least the three brightest AGNs, NGC 1068 (Kinkhabwala 
et al. 2002; Brinkman et al. 2002; Ogle et al. 2003), Circinus 
(Sambruna et al. 2001b) and Mkn 3 (Sako et al. 2000; Bianchi 
et al. 2005b; Pounds & Page 2005) have revealed that their soft 
X-ray emission comes from regions extended on hundreds of 
pc, and the dimension and morphology of this emission closely 
resembles that of the Narrow Line Region (NLR) as mapped 
by the [OIII]/15007 emission line (see e.g. Bianchi, Guainazzi 
& Chiaberge 2006). 

Systematic studies on large samples of FRII Narrow Line 
Radio Galaxies (NLRG) show that the High Excitation Radio 
Galaxies (HEG; Jackson & Rawlings 1997) have a heavily 
obscured continuum (intrinsic column density ~10 23 ~ 24 cur 2 ) 
related to the bright accretion disk obscured by an edge-on 
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oriented thick torus, and a soft unabsorbed tail (Belsole et al. 
2006; Hardcastle et al. 2006; Evans et al. 2006). 
Soft X-ray emission above the extrapolation of the obscured 
component is often detected. The origin of the soft X-ray 
emission in radio galaxies is still matter of debate. Radio-loud 
AGNs are potentially more complex than their radio-quiet 
counterparts, because the jet could also contribute to the X-ray 
emission. 

It was suggested that the soft excess is related to the jet 
emerging from the edge of the dusty torus (Hardcastle & 
Won-all 1999; Belsole et al. 2006). Recently the discovery of 
emission lines produced by the photoionized gas in the soft 
X-ray spectrum of 3C 445 (Grandi et al. 2007; Sambruna et al. 
2007) and 3C 234 (Piconcelli et al. 2008) has opened a new 
and alternative way to interpret the soft excess. 
3C 33 is the right target to investigate the nature of the soft 
excess. It is a genuine type II radio galaxy, optically classified 
as High Excitation Radio Galaxy (HEG) and one of nearest 
HEG in the 3CRR catalogue (Laing Ri ley & Longair 1983) 
located at z=0.0597 jPopescu et al. 19961) . 
Its radio morphology is typical of FRII radio galaxies 
(L 17 g=7.5xl0 25 W Hz 1 sr _1 ), showing simmetric lobes with 
hot spots on kpc scale (Leahy & Perley 1991), and tw o 
simmetric jets in the VLBI image dGiovannini et al. 20051) . 
The southern jet is stronger by a factor 2 at 5 mas. Expressing 
the jet/counterjet ratio (J) in terms of the bulk velocity of the 
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jet (J3) and the angle between the jet and the line of sight (6), Q 
an inclination angle of 75° - 80° can be deduced. Its spatial 
orientation is then favorable to reveal soft X-ray emission 
lines, because a strong de-amplification of the jet radiation is 
expected. 

Moreover 3C 33 is very luminous in the [OIII] band 
(L[o///] ~1.09xl0 42 erg s -1 ) as it exhibits the highest [OIII] 
flux among the nearby NLRGs (z<0.3) of the 3CRR sample. 
Previous nuclear observations of 3C 33 in the X-ray band go 
back to EINSTEIN observatory (Fabbiano et al. 1984) and 
ROSAT (Hardcastle Birkinshaw & Worrall 1998; Hardcastle 
& Worrall 1999) X-ray telescopes. EINSTEIN only provided 
an upper limit to the 0.5-3 keV luminosity, while ROSAT/HRI 
could detect the source in the soft X-ray band, although the 
HRI instrument on-board it could not resolve the source. 
Evans et al. (2006) analyzed the XMM-Newton EPIC PN data 
of 3C 33 showing for the first time the Seyfert 2 nature of the 
hard component. They modeled the soft component with an 
unabsorbed power law. 

Kraft et al. (2007), taking advantage of the unprecedented 
imaging of Chandra, could perform a very detailed study 
of 3C 33. They could resolve the north and south hot spots 
regions. They also studied the nuclear component with a basic 
model consisting of an absorbed power law plus a Fe Ka line 
and a thermal component. Large residuals between 1 and 3 
keV could be explained either by assuming the presence of a 
reflection component or of a less absorbed power law, that they 
interpreted as jet emission on parsec scale. Kraft et al. (2007) 
used thermal gas to model the soft X-ray emission on the base 
of the Chandra imaging that reveals in the soft band (0.3-0.7 
keV) a clearly extended structure around the unresolved hard 
nucleus. A similar elongation was observed in the WFPC1 and 
WFPC2 Hubble Space Telescope imaging, and was interpreted 
by the authors as a tracer of stellar distribution. 
In this work we reanalyzed all the Chandra and XMM-Newton 
observations of the radio galaxy 3C 33 available in the public 
archives, and propose a different interpretation of the nature of 
the soft X-ray emission of 3C 33. 



2. Observations and data reduction 

Chandra ACIS-S observed 3C 33 twice [-20 ks on 2005 
November 8 (ObsID 6190) and 2005 November 12 (ObsID 
7200)]. Spectra and instrument responses were generated using 
the Chandra Interactive Analysis of Observations (CIAO) 3.4 
and the Chandra Calibration Database (CALDB) 3.4.3. The 
source spectrum was extracted from a circular region with 
radius 1.5" centered at the peak of the X-ray emission, while 
background counts were extracted from a 4" radius circle 
in the same CCD, but outside of the extended emission of 
the source. The Chandra count rate from the nucleus in both 
observations is ~0.06 cts s _1 , implying a negligible pileup 
effect. 

In the hard band (2-8 keV) the X-ray morphology of the 
3C 33 nucleus is point-like. An extended structure around the 
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nucleus becomes evident in the soft band (0.5-2 keV), looking 
approximately elliptical, with a major axis of ~ 5", elongated 
in the northeast-southwest (NE-SW) axis, in agreement with 
the analysis of Kraft et al. (2007). 

XMM-Newton observed 3C 33 twice [on 2004 January 4 and 
January 21 (ObsID 0203280101; 0203280301)]. We checked 
both observations but we analyzed only the second pointing (~ 
9 ks) because the first one is affected by high flaring activity. 
This prevented us from obtaining useful information from the 
RGS analysis. 

The EPIC cameras were operated in 'prime full window' 
mode. The data were reduced using the SAS (v. 7.1) with 
standard procedures. 

The light curve over 10 keV was extracted to check high back- 
ground periods. We excluded the time intervals characterized 
by count rate higher than 1.0 cts/s in the EPIC-PN, and 0.35 
cts/s in the EPIC-MOS. After this data cleaning we obtained a 
net exposure of 6.3 ks for the PN and 8.5 ks for the two MOS 
cameras. 

The source and the background spectra were extracted from 
circular regions of 35" radius. The background was chosen 
in the same CCD of 3C 33, but in a region not contaminated 
by the source. The source counts are ~1500 and ~1000 in the 
PN and the two combined MOS, respectively. The response 
matrices were created using the SAS commands RMFGEN 
and ARFGEN. 

Events outside the 0.5-10 keV band were discarded in both PN 
and MOS spectra (Kirsch 2006). 

All spectral fits presented in this paper include absorption 
due to a line-of-sight Galactic column density of N#=3.06 
xl0 20 cirr 2 (Kalberla et al. 2005). 

Errors are quoted at the 90% confidence level (Ay 2 =2.7) for 
one interesting parameter. 

The cosmological parameters used throughout the paper are 
H Q = 70 km" 1 s" 1 Mpc" 1 , Q„, = 0.3, Q A = 0.7 (Spergel et 
al. 2007). For 3C 33, z=0.0597 corresponds to a luminosity 
distance of 263.9 Mpc, and 1" corresponds to a physical scale 
of 1.14kpc. 

The spectral analysis obtained by both telescopes were per- 
formed using the XSPEC 1 1.3 package. 



3. Spectral analysis 

We started our study using Chandra data. The two observations 
were simultaneously analyzed using the spectra extracted from 
1.5" regions around the nucleus. We also considered a larger 
extraction radius of the nuclear source (r=4.5") in order to 
include the extended soft X-ray emission. As expected the 
statistical quality of the spectrum decreased because of the 
higher background. Indeed, very few are the source counts 
in the (0.5-2) keV band from an elliptical region of 5.5" 
semi-major axis and 4" semi-minor axis (excluding the 
nucleus), less than 3% of the total counts in the 1.5" region. In 
physical terms most of the soft-X radiation comes from ~ 2 
kpc around the source. We then decided to consider the smaller 
nuclear region in order to maximize the spectrum signal to 
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noise ratio (S/N). 



3. 1. Hard X-ray spectrum 

The broad band spectrum appeared quite complex. We initially 
decided to fit only the energy range between 3-10 keV in 
order to focus on the hard X-ray continuum excluding the 
soft emission. At first we applied a photoelectrically absorbed 
power-law plus a Gaussian line profile. The restricted energy 
band and the limited statistics do not allow to constrain very 
well the spectral index (the column density and the spectral 
index being strongly correlated), so we decided to fix T=1.7 
following Kraft et al. (2007). This choice was confirmed to 
be the right one by the XMM-Newton analysis (see below). 
As the x L = 102(80) was not completely satisfactory we tested 
for the presence of a reflection component, as in obscured 
AGNs is often required [see for example, Mkn3 (Pounds & 
Page 2005), and 3C445 (Grand i et al. 2007; Sambruna et al . 
2007)]. When a (PEXRAV) dMagdziarz & Zdziarski 1995b 
model was applied the x 2 decreased to 82 for 79 d.o.f., giving 
a better representation of the data (Pf 99.99%) . The reflection 
(R=Q/27r) was let free to vary, while the high-energy cut-off of 
the incident power law was set at 160 keV (Grandi, Malaguti 
& Fiocchi 2006), and the elemental abundances were assumed 
to be solar. We repeated the same excercise leaving F free to 
vary. Also in this case the addition of a reflection component 
improves the fit, even if the strong correlation between the 
reflection component, F and Nw prevent us to substantially 
constrain these three variables. 

The FeKff emission line is not resolved. The rest-frame energy 
is E^ ff =6.38^g gg keV, the upper limit of the intrinsic width is 
cr < 180 eV and the equivalent width is EW & =172+^ eV. 



3.2. Soft X-ray spectrum 

Once a good fit of the hard component was obtained we ex- 
tended our study to the soft part (0.5-3 keV) of the spectrum. 
The simple extrapolation of the hard model introduced a strong 
soft component. Currently, in the literature, at least three 
different interpretations of the soft excess in obscured radio 
loud AGNs are proposed: i) thermal radiation from collisional 
gas (Kraft et al. 2007; Evans et al. 2006); ii) non-thermal 
emission from a jet (Belsole et al. 2006); iii) nuclear scattered 
radiation plus emission lines produced by a circumnuclear 
photoionized gas (Bianchi & Guainazzi 2007; Grandi et al. 
2007; Sambruna et al. 2007; Piconcelli et al. 2008). All these 
three possibilities were tested. 

At first we tried to reproduce it with a thermal model (mekal) 
absorbed by the Galactic column density. This implicitly 
assumes that the soft excess is produced by a collisional 
plasma (with solar abundances). The gas temperature required 
by the fit is kT ~ 0.7, and the model partially reproduced 
the soft excess, ^ 2 (d.o.f.)=138(104), but the inspection of the 
residuals still showed excesses below 2 keV. A second mekal 
was added to verify whether a multi-temperature gas is present, 



but the second temperature shifted to a very improbable value 
larger than 10 keV, and excesses were still present in the soft 
X-ray band (see Fig.l left panel). Note that if the abundances 
are free to vary, they become extremely low, miming a power- 
law. Indeed replacing the "second" thermal emission model 
with a simple power-law (F s free to vary), that represents a 
non-thermal emission, the x 2 improved ^ 2 (d.o.f.)=124(103), 
and the thermal emission is no longer required anymore (the 
mekal model is consistent with zero). Moreover, although this 
was the best fit obtained so far, the model could not completely 
reproduce the soft X-ray data (Fig. 1 right panel). 
Finally the third hypothesis was verified. A soft excess was 
fitted with a power law plus lines, as expected in case of 
circumnuclear photoionized gas. The soft power law slope 
(F s ) was linked to the hard one (F/,), in order to represent 
the scattered nuclear continnum. The fit (see Fig|2]Ze/f panel) 
was very satisfactory ^ 2 (d.o.f.)=103(99). A summary of the 
best-fitting parameters for the hard X-ray spectra is shown in 
Tab[U while the detected lines and relative significance levels 
are listed in Tab|2] 

In order to find indipendent confirmation of this result, we 
explored the XMM-Newton EPIC spectra. The PN and MOS 
results were reassuring. All the EPIC spectra (see Fig|2] right 
panel) were very well fitted with the same best-fit model of the 
Chandra continuum. Interesting enough, the PN data needed 
soft lines at the same energies measured by Chandra. 
The PN best-fit parameters are listed in Tab. [3] The relative x 2 
is 37.5 for 37 degrees of freedom. 



4. Results 

Our analysis displays typical features produced by a cold 
circumnuclear gas: an unresolved (cr < 180 eV) narrow iron 
line (EW ^ 200 eV), and a reflection component (R ^ 2) 
(Risaliti 2002), showing the Seyfert 2 nature of the hard 
spectrum of 3C 33. Both these values are consistent with the 
Chandra results of Kraft et al. 2007 (their Model 1). 
The reflection component was obtained assuming an arbitrarily 
inclination angle of 6 - 60° (PEXRAV parameter cos6> 0.45). 
However even considering a larger inclination angle, as sug- 
gested by the radio data, the value of R does not significantly 
change. This is not surprising as the limited statistics and 
the large number of parameters required by the fit. Although 
the reflection component is slightly higher than that expected 
from an homogeneous dusty torus, it is not unusual and 
often observed in radio-quiet obscured objects. This probably 
indicates the presence of a more complex inhomogeneous and 
structured circumnuclear absorber (Risaliti et al. 2002). 

The spectral analogy between Seyfert 2s and obscured ra- 
dio galaxies seems to be maintained at soft X-ray energies, 
where the best fit requires emission lines overimposed on a 
weak continuum. The lines in Tab. |2] and @] are likely produced 
by H- and He-like ions of light elements such as Oxygen and 
Neon, often detected in obscured AGNs showing strong evi- 
dence of photoionized gas (Bianchi & Guainazzi 2007; Grandi 
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Table 1. Chandra ACIS best-fit parameters. Normalizations n 
are expressed in cirT 2 s _1 keV -1 ; the normalization of the iron 
line is expressed in photons crrr 2 s . 



Table 3. XMM-Newton EPIC best-fit parameters. 
Normalizations n are expressed in cirT 2 s -1 keV -1 ; the 
normalization of the iron line is expressed in photons cirT 2 s _1 . 



r 

R 

Nq 
N« 
n 

"Wiobscured 

E Fe (keV) 
F Fe 

EW (eV) 




1.7* 

9 i q+0.48 

(3.06X10 20 ) cm" 2 
(5.28^'°J) x 10 23 cm" 2 
(1.83:0^) x 10- 3 
(4.48+ 3 - 43 ) x 10- 6 

t OO+0.05 

"•^-O.oe 
(2.491| n)x 10- 5 
172 +84 

1 ' -76 


r 

R 

N G 
N h 
n 

^unobscured 

E Fe (keV) 

F Fe 

EW (eV) 




1 69 +018 

9 7+O.8I 
A- '-0.60 

(3.06xl0 20 ) cm" 2 
(7.14 + | ^) x 10 23 cm" 2 
(1.83*$) x 10- 3 
(8.1 +4 -° 8 ) x 10- 6 
6 40 +013 

"•^"-0.08 

(3.65^ 3 ) x 10- 5 
206i|* 


^(d.o.f) 




103(99) 


r(d.o.f) 




37.5(37) 


Unabsorbed luminosity and absorbed flux 

Li-w fcv(ergs-') 7.31 XlO 43 
Flux2-m kev (erg s~ ] cm -2 ) 2.67 xlO" 12 


Unabsorbed luminosity and absorbed flux 

£2-10 fevCergs- 1 ) ' 7.53X10 43 
Flux 2 -to kev (erg s - ' cnT 2 ) 2.64 xlO -12 



(*) value fixed 

Table 2. Emission lines detected in the Chandra ACIS spec- 
trum. E c is the centroid line energy (keV) in the source rest 
frame; F is the line flux in units of 10~ 6 photons crrr 2 s _1 ; Pf 
is the significance of the fit improvement after including the 
Gaussian line in the model, based on the F-test. 



Energy 


Intensity 


Pf 


Tentative IDs 


76 +0 02 

"■'"-0.13 

n cn+o.o4 

1 30 +004 
L -002 


9 QQ + 1.16 

Z - 7J -1.02 
9 c 1 +0.85 
^-^ -1.20 

1.41±0.46 


> 98.9% 
> 99% 

> 99.97% 


OVII RRC, Fe XVII 3s-2p 
Ne IX rif 
Mg XI Heo- 



Table 4. Emission lines detected in the XMM-Newton EPIC 
spectrum. E c is the centroid line energy (keV) in the source 
rest frame; F is the line flux in units of lCT 6 photons crrr 2 s -1 ; 
Pf is the significance of the fit improvement after including the 
Gaussian line in the model, based on the F-test. 



Energy 


Intensity 


Pf 


Tentative IDs 


0.71+0.03 

93+ 03 

"•"-0.05 

1 26 +004 

i - z '"-0.17 


7.89! 2 « 
2.35±2.2 


> 99.8% 

> 99% 

> 88% 


OVII RRC, Fe XVII 3s-2p 
Ne IX rif 
Mg XI Heo- 



et al. 2007; Sambruna et al. 2007; Piconcelli et al. 2008). 
Although the moderate ACIS energy resolution does not allow 
to draw strong conclusions, we stress that the centroid of the 
lines occurs where two strong features produced by photoion- 
ized gas are expected, i.e. the OVII recombination feature and 
the NelX triplet lines. Incidentally we note that the intrinsic 
nuclear luminosity of 3C 33 (Li_iooond ~2xl0 44 erg sec -1 ) is 
comparable to those observed in typical Seyfert 2s, for which 
the higher quality of the data has permitted detailed spectral 
studies of the soft X-ray diagnostics (see, for example, Pounds 
& Page 2005, Sako et al. 2000ab). Approximate computations 
indicate that the energy budget of the AGN is sufficient to pro- 
duce a range of values of the ionization parameter (£ = L '-™°'-" f ) 
consistent with the photoionized scenario (Kinkhabwala et al. 
2002). If the soft X-ray emission is related to the Narrow Line 
Region (as in Seyfert 2s), assuming a minimum (R m ,„ ~1 pc) 
and a maximum (R„ MA ~2 kpc, our nuclear extraction region) 
distance of the gas, a typical NLR electron density n e ~ 10 3 
cm- 3 implies log£ spanning between -2 and 4 erg cm s _I . 

5. The origin of the soft X-ray emission 

5. 1. Imaging 

X-ray spectroscopy provides evidence for the presence of soft 

X-ray emission lines in the radio galaxy 3C 33. 

In spite of the limited spectral power of the CCDs we could 
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Fig. 1. Residuals corresponding to the two different soft excess 
models: left two mekal and right a scattered power-law. 
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Fig. 2. Chandra spectra (left panel) and XMM spectra (right 
panel) of 3C 33. PN, MOS 1 and MOS 2 spectra are 
represented with black, red and green crosses, respectively. 
Residuals are in terms of standard deviations of the observed 
data from the model. 

detect features that are plausibly emitted by a photoionized gas 
(see Tab. and Tab. @). 

The analogy with Seyfert 2 galaxies is straightforward. Narrow 
Line Radio Galaxies would share similar continua and environ- 
ments with their radio-quiet counterparts. In Seyfert 2 galaxies 
the coincidence in morphology and dimension of the soft X-ray 
and the [OIII]/i5007 emission suggests the Narrow Emission 
Line Region as the location of such gas. It is then natural to 
search for a similar X-ray and optical spatial correlation in 
3C33. 

For this reason the narrow-band HST ACS/WFC1 optical im- 
age, centered at redshifted [OIII]/i5007 (filter FR551N), was 
retrieved from Multi Mission Archive at STScI (MAST) and 
was aligned with the Chandra image as described in Bianchi, 
Guainazzi & Chiaberge (2006). 

Fig.[3]shows the Chandra soft X-ray contours (< 2 keV) super- 
imposed on the HST [OIII]. The coincidence in extension and 
morphology between the two emissions is impressive. 

5.2. Further hints from optical spectroscopy 

Further support in favor of the photoionization scenario comes 
from optical spectroscopy. 

Several line diagnostic diagrams have been proposed to dis- 
criminate between collisionally ionized plasma (e.g. starburst 
emission) and photoionized plasma irradiated by a central en- 
gine (see for example Miller et al. 2003). One of these, based 




Fig. 3. Superposition of the Chandra soft X-ray (<2 keV) con- 
tours on an HST image taken through a linear ramp filter at 
redshifted [OIII]/i5007. The sign "x" indicates the centre of the 
hard X-ray source, north is up, east to the left. The X-ray im- 
age was smoothed with a Gaussian of FWHM ~ 6 pixels. The 
contours correspond to four logarithmic intervals in the range 
1-60% of the peak flux. 

on the [Om]/Hp and [Nil] /Ha ratios (Fig. 3 of Miller et 
al. 2003) could be used to test the nature of the gas emis- 
sion of 3C 33. Using the nuclear reddening-corrected emis- 
sion line luminosities reported by Koski (1978) we estimated 
log([OIII]/Hj0)~ 1.1 and log([NII]/Har)~ -0.17. Such values 
place 3C 33 in the region populated by the photoionized AGNs. 
Baum et al. (1992) performed an accurate study of the ex- 
tended kpc emission-line region (EELR). They also found a 
large value of [OIII] /H/3 ratio that they interpreted as produced 
by photoionized gas with a high ionization parameter. 
This result strengthens the early conclusion of Baum & 
Heckman (1989) who demonstrated that the number of photons 
needed to photoionize the emission-line nebula of this source 
equals the number of ionizing photons coming from the central 
engine. 

In addition, on the base of the kinematic classification of 
powerful radio galaxies proposed by Baum et al. (1992), 
3C 33 belongs to the rotator class. Rotators, characterized 
by a systematic gradient in the velocity of the emission-line 
gas, are nearly always associated with powerful radio sources 
(Lradio >10 42 erg s _1 ) having luminous high-ionization nuclear 
emission-line region. 

The kinematic study of the gas also shows that the gas rota- 
tion axis and the radio axis are aligned within 10° (Baum et al. 
1992). Moreover the rotation axis of the gas is not aligned with 
the stellar rotation axis (Heckman et al. 1985), disfavoring the 
hypothesis of a collisional (starburst) gas emission. 

6. Conclusions 

An X-ray analysis of the Narrow Line Radio Galaxy 3C 33 
is presented. All the Chandra and XMM-Newton observations 
available in the public archives are used. Our study confirms 
the Seyfert 2 nature of this source. In agreement with previous 
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works we confirm that the broad-band spectrum is complex and 
largely dominated by strong obscuration (N H ~10 23 cm -2 ). 
A prominent Fe Ka line at 6.4 keV is also present. Its EW of ~ 
200 eV is in agreement with being produced by optically-thick 
circumnuclear matter. 

We advocate photoionized gas for an interpretation of the soft 
X-ray emission. The clues favoring this hypothesis are: 

1 . both Chandra and XMM-Newton soft X-ray spectra can be 
better fit with highly ionized lines and unobscured scattered 
power-law, rather than with thermal components; 

2. a rough estimate of the ionization parameter expected in 
the case that the emitting line gas is located in the NLR, is 
consistent with a photoionized scenario; 

3. there is a remarkable spatial coincidence between the 
extended soft X-ray and the [OIII]/15007 emission; 

X-ray findings are also consistent with pieces of evidence 
coming from other wavelengths such as: 

4. the nuclear line ratios [OIII]/H/3 and [Nil] /Ha put 3C 33 
in the right position of the line diagnostic plane of Miller et 
al. (2003). It is among the AGNs, far away from the region 
populated by starburst galaxies; 

5. the extended kpc emission region has a line rich optical 
spectrum with large [OIII]/Hfi ratio, typical of photoion- 
ized gas; 

6. the estimated photons required to photoionize the emission- 
line nebula (in the optical band) are provided entirely by the 
nuclear engine; 

7. the kinematical properties of the extended gas of 3C 33, 
based on optical spectroscopy, are typical of powerful radio 
galaxies with high-ionization nuclear emission; 

8. the gas rotation axis and the radio axis are aligned within 
10°, but the former is not aligned with the stellar rotation 
axis. 
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